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Structures of excited states in are investigated with a method of j3-'y constraint antisym- 
metrized molecular dynamics in combination with the generator coordinate method. Various excited 
states with the developed 3a-cluster core structures are suggested in positive- and negative-parity 
states. In the positive-parity states, triaxial deformed and linear-chain structures are found to con- 
struct excited bands. Interestingly, ^"Be-l-a correlation is found in the cluster states above the 
^''Be-l-a threshold energy. 
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^ . I. INTRODUCTION 

■ Owing to the experimental progress of unstable nuclei, various exotic structures have been discovered in neutron- 
d ] rich nuclei where excess neutrons play important roles. In light neutron-rich nuclei, cluster structure with excess 
neutrons is one of the hot subjects in experimental and theoretical studies, where searching for new cluster states in 
excited states has been performed. 
0^ , In this point of view, neutron-rich C isotopes such as ^^C are an interesting subject because they have excess 
neutrons compared with ^^C, which is already well known to show various structures owing to the coexistence of 
(~| \ shell-model and cluster features. The ground state of ^^C has mainly a shell- model feature of the p3/2-subshell 
closed configuration, whereas, in the excited states above the 3a threshold energy, various configurations of 3a- 
cluster structures were suggested in many theoretical works For example, the 3^ state has been discussed 

in association with an equilateral-triangular structure of three a clusters. The a condensation of weakly interacting 
^ ] three a clusters suggested in the and 2 J states is another cluster aspect attracting great interest recently 
—I, Moreover, a linear-chainlike structure with an obtuse-angle-triangular 3a configuration was suggested in the Og state 
(qI. ITo|. Therefore, it is expected that ^^C with two excess neutrons may also exhibit rich phenomena in the excited 
states. In particular, from an analogy of "'^^C, various cluster structures may appear in ^^C as well as shell-model 
structures. 
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' For the excited states of ^^C, there are many experimental studies that indicate the appearance of cluster states 



^\ ■ such as ^"Be-l-a cluster states [i3 - [Tq |. Also on the theoretical side, cluster structures of ^'^C were suggested 17 

■ In the study of ^"^C with 3a+2ri cluster models, Itagaki et al. predicted that an equilateral-triangular structure of 
. the well-developed three a clusters surrounded by excess neutrons is formed constructing K"^ = 0+ and K"^ = 3~ 

■ rotational bands in the excited states. They argued that these two bands originate in a rather rigid 3a structure 
, stabilized by the excess neutrons in the molecular orbitals. This mechanism can be interpreted as a realization of 
• the a crystallization in dilute nuclear medium. Another interesting problem is whether a linear-chain 3a structure 

^ I with two excess neutrons exists in the excited states of ^^C. Although there were discussions on the possibility of the 
' ^ • linear-chain structure in neutron-rich C isotopes fl^ [isl. [l9| . there is no clear conclusion for stability of the linear-chain 
r^J ] structure in ^"^C. Thus, the excited states of ^^C are attracting much interest recently, and therefore, systematic study 
of the ground and excited states is required. 

Our aim is to investigate the ground and excited states of ^^C while focusing on cluster features. To clarify the role 
of excess neutrons in neutron-rich nuclei, it is helpful to consider analogies and differences of cluster features between 
^'^C and ^^C. It is also a challenging problem to search for new types of structures in ^^C. In such a systematic 
study, it is rather important to apply a framework which is free from model assumptions of clusters because the 
formation of three a clusters is not obvious in ^"^C. Therefore, we adopted a method of antisymmetrized molecular 
dynamics (AMD), which is a framework without any assumption of the existence of clusters. In the present study, 
we apply a combination of the /3-7 constraint AMD and the generator coordinate method (GCM), which we call the 
/3-7 constraint AMD -I- GCM. This method has already been proved to be a powerful approach to describe various 
structures such as cluster and shell- model-like structures In particular, it is useful for the systematic study of 
cluster states in excited states owing to the superposition of basis AMD wave functions on the two-dimensional /3-7 
plane. Therefore, it is suitable for the study of ^""^C. 

This article is organized as follows. In Sec. |lll we explain the framework of the f3-j constraint AMD -I- GCM. The 
calculated results are shown in Sec. IIIII In Sec. IIVI we discuss the effect of excess neutrons and compare the present 
results with earlier works. Finally, in Sec. [Vj a summary and an outlook are given. 
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II. FRAMEWORK OF THE /3-7 CONSTRAINT AMD + GCM 

The frameworks of AMD are described in detail, for example, in Refs. [2T| - [23 |. In this article, we adopt a version 
of AMD, the 13- j constraint AMD [2^, in which wc perform the variation with the constraint on the quadrupole 
deformation parameters, /3 and 7. 



A. Wave function of AMD 



In the method of AMD, a basis wave function of an A-nucleon system |$) is described by a Slater determinant of 
single-particle wave functions as 



1$) 
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det{|v3i),--- ,\ipA)}- 



(1) 



The z-th single-particle wave function \ipi) consists of the spatial part 10,), spin part and isospin part \Ti) as 

\^,) = mxr)\n). (2) 
The spatial part \(f)i) is given by a Gaussian wave packet whose center is located at Zi/^/v as 



exp 



(3) 



where v is the width parameter and is taken to be a common value for all the single-particle Gaussian wave functions 
in the present work. The spin orientation is given by the parameter ^i, while the isospin part \Ti) is fixed to be up 
(proton) or down (neutron). 



\X^)=^^t\t)+^^l\i), 

In) = \p) or |7i). 



(4) 
(5) 



In a basis wave function |$), {X} = {Z , ^} = {Zi, ^1, Z2, ^2, • • • , Za, ^a} are complex variational parameters and 
they are determined by the energy optimization. 



B. Parity and angular momentum projections 

We project the AMD wave function onto parity and angular momentum eigenstates by using the parity projection 
operator and the angular- momentum projection operator Pj^jj^- The parity projection operator P^ is defined as 

. (6) 

where P is the parity operator. The angular-momentum projection operator P^jj^ is defined as 



Pmk ^ J dnDi*^.{n)Rm, (7) 

where = (a,/3,7) are the Euler angles, D'ljj^{^) is the Wigner's D function, and R{Vl) = ^-^^j^ ^-ipjy ^-i-fj^ jg ^j^g 
rotation operator. 

We perform the variation for the parity projected wave function 1$*) defined as 

|$±)=P±|$). (8) 

After the variation, we project the obtained wave function onto the total-angular-momentum eigenstates. It means 
that the parity projection is performed before the variation, and the total-angular-momentum projection is carried 
after the variation. 
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C. Constraint variation on quadrupole deformation 

To describe various cluster and shell- model structures which may appear in the ground and excited states of ^^C, 
we perform the energy variation with the constraints on the quadrupole deformation parameters, /3 and 7. 
The deformation parameters, /? and 7, arc defined as 

^,i„..;f<^2^. (10, 

s I + m + {€-)) . (11) 

Here, (O) represents the expectation value of the operator O for an intrinsic wave function |$). x, y, and z are 
the inertia principal axes that are chosen as (y^) < (i^) < (i^) and (xy) = (yz) ~ (zx) ~ 0. To satisfy the latter 
condition, wc also impose the constraints {xy) j E? = {yz) j ~ {zx)/B?' = 0. To obtain the energy minimum state 
under the constraint condition, wc add the constraint potential T^onst to the total energy of the system in the energy 
variation. The constraint potential K;onst is given as 



Konst =m [(/? COS 7 - /3o COS 70)^ + (/3 siu 7 - /3o sin 70 )^ 



+1I2 



Here, ?yi and 772 take sufficiently large values. That is, we minimize the constrained energy i?^nst defined as 



(12) 



-^const = ^(j)± |([)±^ ^ Kionstj (13) 

where H is the Hamiltonian, to determine {X}. After the variation with the constraints, we obtain the optimized 
wave functions |<i>^(/3o, 7o)) for each set of parameters, (/3,7) = (/3o,7o)- 



D. Generator coordinate method 



In the calculations of energy levels, we superpose the parity and total-angular-momentum projected AMD wave 
functions P/^j^|<i>^(/3, 7)). Thus, the final wave function for the state is given by a linear combination of the basis 
wave functions as 

lO = EE/"(/'-7.,^)Pi^i.i*^(A,70)- (14) 

K i 

The coefficients fn{(3i,'yi, K) are determined using the Hill- Wheeler equation 

5 [{^i^imi^) = 0. (15) 

This means the superposition of multiconfigurations described by parity and total-angular-momentum projected AMD 
wave functions. In the limit of sufficient basis wave functions on the /3-7 plane, it corresponds to the GCM with the 
two-dimensional generator coordinates of the quadrupole deformation parameters, /3 and 7. 

E. Hamiltonian and parameters 

The Hamiltonian H consists of the kinetic term and effective two-body interactions as 

H = Y.U~fG + Y, + E + E , (16) 

i i<j i<j i<j 
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pin-orbit ^ „ j WCoulomb 

the central force, we use the Volkov No. 2 interaction [2 



where V'j^''"*'''^', y^^P"'""' ^ and ^^ouiomb g^.^ ^j^g central force, spin-orbit force, and Coulomb force, respectively. As 



V, 



central 



E 



Vk exp 



{W + BP, - HPr - MP^Pr), 



(17) 



where wi = —60.65 MeV, V2 = 61.14 MeV, ai = 1.80 fm, and 02 = 1.01 fm. For the spin-orbit part, we used the 
spin-orbit term of the GSRS interaction [26|, which is a two-range Gaussian with a projection operator P('^O) onto 
the triplet odd state. 



spin-orbit 



E 

k=l 



Uk exp 



pCo)l-s, 



pCo) 



1+P,l + Pr 



(18) 
(19) 



where bi = 0.600 fm and 62 = 0.477 fm. Here, P, and Pr are the spin and isospin exchange operators, respectively. 

We take the same interaction parameters as those in Ref. [2^, i.e., the Majorana exchange parameter M = 0.6 
{W = 0.4), the Bartlett exchange parameter B = 0.125, and the Heisenberg exchange parameter H = 0.125 in the 
central force, and ui = —1600 MeV and U2 = 1600 MeV in the spin-orbit force. All these parameters are the same as 
those adopted in the studies for ^Be [13, ^°Be [lil, and C isotopes, such as ^^C, "C, and ^^G [IT 
small modification in the strength of the spin-orbit force. They are adjusted to reproduce the a + 
W ~ I — M , ai, 02), binding energy of the deuteron (B ~ H), and a + n phase shift (ui = —U2, 



the slightly weaker strengths of the spin-orbit force than — ui = U2 = 2000 MeV adopted in Refs 
fit the energy of ^^C 0. 

For the width parameter of single-particle Gaussian wave packets in Eq. ([3]), we used the value 
which is determined from a variational calculation for the ground state of ^Be in Ref. 
as those in the studies for Be isotopes and G isotopes [13, [H, [l^l ■ 



], [18[, except for a 
a phase shift (M, 
6], 62L We adopt 

[B, [H m m to 



0.235 fm" 



27| . This value is also the same 



III. RESULTS 



We applied the /3-j constraint AMD + GGM to ^ G. In this section, we show the calculated results. 



A. Energy surfaces 



We performed variational calculations with the /3-7 constraint at 196 mesh points of the triangle lattice on the /3-7 
plane. Energy surfaces as functions of /3 and 7 are obtained. The calculated energy surfaces on the /3-7 plane for 
positive-parity states and negative-parity states are shown in Figs. [T] and [21 respectively. 

In Fig.[TJ the top panel shows the energy of the positive-parity states before the total-angular- momentum projection, 
and the bottom panel shows the results for the 0"*" states calculated by the total-angular-momentum projection after 
the variation. We call the former the positive-parity energy surface and the latter the 0"*" energy surface. The minimum 
point of the positive-parity energy surface is at (/3 cos 7, /3 sin 7) = (0.00, 0.00), which indicates a spherical shape. After 
the total-angular- momentum projection, the minimum point of the 0^ energy surface becomes (^ cos 7,/? sin 7) = 
(0.23,0.04). It means that the deformation of the energy minimum state changes from the spherical shape before 
the total-angular-momentum projection to the triaxial shape after the projection, because higher correlations beyond 
mean field are incorporated by the total-angular-momentum projection. In a largely deformed region, there exists a 
flat region around {(Hcosj, fSsm-f) = (0.80 — 0.90, 0.00) in the positive-parity energy surface. After the total-angular- 
momentum projection, a local minimum at (/3cos7,/?sin7) ~ (0.90,0.04) emerges from this flat region. As we show 
later, a rotational band with the large prolate deformation is constructed by wave functions in this region after the 
GGM calculation. 

The negative-parity energy surface and the 3~ energy surface are displayed in the top and bottom panels of Fig. [51 
respectively. The minimum point of the negative-parity energy surface is at (/3 cos 7,/? sin 7) = (0.08,0.04), and that 
of the 3~ energy surface is at (/3cos7,/3sin7) = (0.20,0.09). Again, the deformation of the energy minimum state 
changes from an almost spherical shape to a triaxial one before and after the total-angular momentum projection. In 
the largely deformed prolate region, the behavior of the energy surface for the negative-parity states is different from 
that of the positive-parity states. Along the 7 = 0° line, the energy rapidly increases and there is no flat region in 
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FIG. 1: Energy surfaces of C on the /3-7 plane. The upper shows the energy for the positive-parity states before the 
total-angular-momentum projection and the lower shows that for the 0^ states after the total-angular-momentum projection. 

the negative-parity surface. Even after the total-angular-momcntum projection, no flat region is found around the 
largely deformed prolate region as seen in the 3~ energy surface. 

B. Structures on the /?-7 plane 
In this section, we describe intrinsic structures obtained by the constraint AMD. 

We analyze the spatial configurations of the Gaussian centers {Zi, Z2, ■ ■ ■ , Zj^} and the distributions of proton 
density pp and neutron density p„ of each intrinsic wave function |$(/3,7)). We also investigate the neutron-proton 
density difference pn — Pp to observe excess neutron behaviors. To demonstrate density distributions, we show the 
density p integrated along the y-axis as 

pix,z)= J dypir), (20) 

p{r) EE ($(/3,7)| ^5(r - f,)|$(/3, 7))- (21) 

i 

First, we discuss the intrinsic structures of positive-parity states of ^"'C. The density distributions of the intrinsic 
wave functions for positive-parity states are illustrated in Fig. [31 The energy minimum state at (/3 cos 7, /3 sin 7) = 
(0.00,0.00) in the positive-parity energy surface shows almost spherical density distributions as seen in Fig. [31[a). In 
this wave function, the centers of the single-particle Gaussian wave packets gather near the origin. That is, this state 
has no spatially developed cluster structure and it is almost equivalent to the shell- model state with the p3/2-subshell 
closed-proton configuration and the p-shcll closed-neutron configuration. The density distribution for the minimum 
point (/3cos7,/3sin7) — (0.23,0.04) in the 0+ energy surface is shown in Fig.|3l^b). This wave function is found to 
be the dominant component of the ground state obtained by the GCM calculation as is shown later. In this state, 
an a cluster core somewhat develops compared with the state (a) for the energy minimum before the total-angular- 
momentum projection. However, single-particle Gaussian wave packets still gather around the origin; therefore, this 
state is regarded as the intermediate between the shell-model structure and the cluster structure. This feature is also 
reflected in the expectation value of squared intrinsic spin of protons (5*^) = 0.77, which is smaller than the ideal 
value 4/3 for the p3/2-subshell closed configuration and larger than the value for a 3a-cluster configuration. 

Figure |3Jc) indicates the density distribution for a triaxial deformed state. In this state, to be clear from the 
proton density and also from the expectation values of the squared intrinsic proton spin (Sp) = 0.15, three a-cluster 



6 



[MeV] 




0.0 0.2 0,4 0.6 0.3 1.D 1.2 1.4 
j8 cos T 



FIG. 2: Energy surfaces of C on the fi-j plane. The upper shows the energy for the negative-parity states before the 
total-angular-momentum projection and the lower shows that for the 3^ states after the total-angular-momentum projection. 

cores develop. Excess neutrons occupy the sd-like orbitals distributing along the bisector of the vertical angle of the 
isosceles-triangle which consists of three a clusters. These neutron orbitals make this state triaxial. After the GCM 
calculation, a = O"*" band and its = 2+ side band dominated by this triaxial state are obtained. 

In the large deformation region, three a-cluster cores develop well in ^^C. Various configurations of three a clusters 
appear, depending on the deformation parameters, /? and 7, as seen in Figs. [31[d), (e), and (f) for typical density 
distributions of oblate, triaxial, and prolate deformed states, respectively. It is found that the equilateral-triangular 
structure, obtusc-anglc-triangular structure, and linear-chain structure of three a clusters arise in the oblate state 
(d), the triaxial state (c), and the prolate state (f), respectively. In these states, the expectation values of squared 
intrinsic proton spin (Sp) are smaller than 0.10, which is almost consistent with the value for three (Os)'^ a cores. 
It is interesting that the excess neutrons distribute around two of the three a clusters in the triaxial state (c) and 
prolate state (f), which indicates ^"Bc+a correlations. In the viewpoint of two-center ^°Be-fa cluster structure, the 
degree of the freedom for rotation of the deformed ^"Be cluster is described on the P-j plane. We mention here that 
the linear-chain structure (f) is the base of the flat energy region at the large prolate deformation. After the GCM 
calculation, the linear-chain state (f) constructs a rotational band. 

Next, we discuss the intrinsic structure of negative-parity states. The density distributions for negative-parity 
states are illustrated in Fig. 21 The density distribution of the intrinsic wave function at the minimum point 
(/3cos7, /3sin7) = (0.08,0.04) in the negative-parity energy surface is shown Fig. IH^a). In this wave function, the 
centers of the single-particle Gaussian wave packets gather near the origin and there is no spatially developed cluster 
structure. In this state, the expectation value of squared intrinsic spin of protons is 1.24, which is close to the ideal 
value 4/3 for the p3/2-subshell closed configuration, while that of neutrons is 0.85, which deviates from the p-shell 
closed configuration value 0. Therefore, this state is interpreted as the shell-model structure having Ip-lh excitations 
on the neutron p shell with the almost p3/2-subshell closed-proton configuration. Compared with an almost spherical 
shape of the state (a), an a cluster core begins to develop slightly at the minimum point (/3 cos 7,/? sin 7) — (0.20,0.09) 
in the 3' energy surface (Fig. UJJb)). However, single-particle Gaussian wave packets still gather around the origin, 
and this state is regarded as the intermediate between the shell-model structure and the cluster structure. 

Figures lll^c), (d), and (e) show intrinsic density distributions for typical deformed states — oblate, triaxial, and 
prolate, respectively. It is found that the isosceles-triangular structure and obtuse-angle-triangular structure of three 
a-cluster cores arise in the oblate state (c) and triaxial state (d), respectively, as shown in the proton densities. These 
structures are similar to those obtained for the positive-parity states. 

However, the structure of the prolate deformed states is somewhat different from the positive-parity states. In the 
prolate state (e), the proton density pp has three peaks, which indicate three a-cluster cores. However, the expectation 



7 



Pv 



pn 



Pn — Pp 






i/im J 




0.625 


1 


0.575 




0.525 




0.475 




0.425 




0.375 




0.325 




0.275 




0.225 




0.175 




0.125 




0.075 


! 


0.025 




-0.025 




-0.075 



FIG. 3: (Color online) Density distributions of the intrinsic wave functions for the positive-parity states of ^^C. The proton 
density pp, neutron density p„, and difference between the neutron and proton densities pn — Pp are illustrated in the left, middle, 
and right columns, respectively. The density distributions of the intrinsic wave functions at (a) (/3 cos 7, /3 sin 7) = (0.00,0.00), 
(b) (^g cos 7, /3 sin 7) = (0.23,0.04), (c) (/3 cos 7, ^ sin 7) = (0.45,0.17), (d) (/3 cos 7, /J sin 7) = (0.25,0.35), (e) (/3 cos 7, /3 sin 7) = 
(0.78,0.22), and (f) (/3 cos 7, /3 sin 7) = (0.93,0.04) on the /3-7 plane are shown. The size of the box is 10 x 10 fm^. 



values of squared intrinsic proton spin (Sp) = 0.33 indicates that the prolate state (e) contains components of the (Os)* 
a breaking. Moreover, three a clusters shows a bending structure in stead of the straight-line 3a configuration seen in 
the linear-chain structure (Fig.[3ljf)) of the positive-parity prolate state. That is, the ^''Be cluster rotates considerably 
though the prolate state (e) has almost the axial symmetric shape caused by the constraint parameters close to the 
7 = 0° axis on the plane. This indicates that the linear-chain structure is less favored in the negative-parity 
states than in the positive-parity states. It is consistent with the features of the (3-"f energy surfaces, where there is 
no plateau in the negative-parity states as mentioned before. 

Thus, on the /3-j plane for the ^^C system, we obtained various structures such as the shell-model structure in the 
small (3 region and the cluster structure in the large /3 region. In particular, various spatial configurations of three a 
clusters such as the linear-chain and the equilateral-triangular structures are obtained as a function of the triaxiality 
7. These basis wave functions are proper for the study of where cluster and shell-model aspects are expected to 
coexist. 



C. Energy levels 

To calculate the energy spectra of "^^C, we superposed the wave functions obtained by the constraint AMD at 
196 mesh points on the plane by using the GCM. 

First, we describe the results for the positive-parity states. We show the calculated positive-parity energy levels in 
Fig. [5]as well as the experimental levels. In the three columns from the left, we display the experimental energy levels 
for all the positive-parity assigned states [2§| . In the fourth column, we show the states observed recently in ^^Be+a 
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FIG. 4: (Color online) Density distributions of the intrinsic wave functions for the negative-parity states of ^^C. The proton 
density pp, neutron density p„, and difference between the neutron and proton densities pn — Pp are illustrated in the left, 
middle, and right columns, respectively. The density distributions of the intrinsic wave functions at (a) (/3 cos 7, /3 sin 7) = 
(0.08,0.04), (b) (/3 cos 7, sin 7) = (0.20,0.09), (c) (/3 cos 7, ^ sin 7) = (0.25,0.35), (d) (/3 cos 7, ^ sin 7) = (0.60,0.17), and (e) 
(/3cos7, /3sin7) = (0.93,0.04) on the 13-j plane are shown. The size of the box is 10 x 10 fm^. 



breakup reactions [T3-[TB|- The lower three of them are the states with the strong population of the ^'^Bc{Oi)+a 
and the upper three are those of the ^^Bc{2^)+a. In the theoretical results, we classified the calculated states in 
five groups by analyzing the components of the basis wave functions and E2 transition strengths. The calculated E2 
transition strengths and root-mean-square radii are listed in Tables U and |TT1 respectively, along with experimental 
data [H,!!!. 

The calculated 0^ and 2^ states constitute the ground band. The calculations reproduce well the experimental 
data for the properties of the ground band such as B{E2, 2^ — Of) and the root-mean-square radius of the 0^ state. 
The results show that the basis wave function at the minimum point (/3 cos 7, /3 sin 7) = (0.23,0.04) in the 0"*" energy 
surface dominates these ground-band states. For instance, the overlap with the 0^ state is 97%. As discussed before, 
this dominant basis wave function has no spatial development of three a clusters but the intermediate structure 
(Fig.lSb)). 

The calculated Oj, 2^, and 4^ states in the = Oj band have large overlap with the AMD base at 
(/3 cos 7, /3 sin 7) = (0.45,0.17) as 78% overlap in the Oj state. The 2^ ,3^ , 4:2, states also contain domi- 

nant components of the same basis wave function. As mentioned before, this basis wave function has the triaxial 
intrinsic structure with three a clusters and excess neutrons occupying the sd-like longitudinal orbitals (Fig. |3tc)). 
Therefore, it is natural to interpret these states as the K'^ = 0+ band and the A'^ = 2+ side band which are con- 
structed from the triaxially deformed intrinsic state. It is an interesting suggestion that the side band is built on the 
O2 band of ^"'C owing to the triaxial shape caused by the excess neutron distributions around the 3a core. The cal- 
culated E2 strengths for the intraband transitions are comparable to that in the ground band. Because there is some 
mixing of K components in these bands, E2 transitions from 4j and 5^ states are scattered. Although the excitation 
energy of the band- head Oj state is much higher than the experimental energy of the Oj state, the calculated Oj, 
2^, and 4^ states may correspond to the experimental 0^, 2 J, and 4^ states when we measure the energies from the 
^''Be-l-a threshold. To conclude the correspondence with the experimental levels, we need more detailed experimental 
data such as transition strengths. 

For the calculated 0^, 2g, and 4^ states, we found that the main component comes from the AMD wave function at 
(/3cos7, /3sin7) = (0.93,0.04) as 64% overlap in the 0^ state. As discussed before, this dominant wave function has 
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FIG. 5: Energy levels of the positive-parity states in ^*C. Four columns on the left are the experimental data and six columns 
on the right are the calculated results. The dotted lines on the left and right show the experimental and theoretical ^'^Be-l-a 
threshold energies, respectively. 



the linear-chain structure of three a-cluster cores with the ^'^Hc+a correlation (Fig. [3Jf)). In the viewpoint of two- 
center ^"Be+a cluster structure, the rotational mode of the deformed ^°Be cluster is at least partially incorporated 
in the GCM on the /3-7 plane. For example, the basis wave function at (/3 cos 7, /3 sin 7) = (0.78,0.22) in the finite 
7 region shows a tilting ^°Be cluster configuration (Fig. [3l^e)). In the results of the GCM calculations, the 0'^ has 
only 9% overlap with this tilting ^°Be configuration, which is much smaller than the overlap with the linear-chain 
structure at (/? cos 7, /3 sin 7) = (0.93, 0.04). In other words, the GCM amplitudes concentrate on the straight-line ^°Be 
configuration rather than the tilting one. As a result of the large deformation of the linear-chain structure of three 
a clusters, the E2 transitions in this band are strong as B{E2,2j 0^) = 67.2 e^fm^ and B{E2,4j 2^) = 88.3 
e^fm"^. Because the ratio B{E2,4^ — )• 2^)/i?(_E2, 2^ —J- 0^) = 1.31 is consistent with that for the rigid rotor model 
10/7, this band is considered to be the rotational band of the linear-chain structure with the ^'^Be-t-a correlation. 
These linear-chain states may be the candidates for the states observed recently in ^°Be-|-Q; breakup reactions. When 
the energies relative to the ^'^Be+a threshold are considered, these experimental states are observed in the energy 
region similar to the theoretical linear-chain band. The significant branching ratios of ^"Be+a decays of these states 
may support the ^"^Be-I-Q! cluster structure. We should comment here that such a rigid rotational band of the 3a 
linear-chain structure does not appear in ^^C. In the stabilization of the linear-chain structure in ^'*C, excess neutrons 
are expected to play an important role. Detailed discussion of its mechanism is given in the next section. 

Other states have no specific structure and are difficult to be classified as band members. The observed 1+ states are 
missing in the present calculations, because spin-aligned configurations are unfavored; therefore such the states cannot 
be described in the constraint variation where only the lowest energy solution is obtained at each deformation. To 
describe the low-lying 1"*" state, it is necessary to introduce other constraints, for example, on the expectation value 
of intrinsic spin. 

Next, we describe the results for the negative-parity states. We show the calculated energy levels of the negative- 
parity states in Fig.[6]as well as the experimental levels. In the four columns from the left, we display the experimental 
energy levels of the negative-parit y as signed states [2^ . In the fifth column, we show the energy levels recently observed 
in ^"Be+a breakup reactions [13l - [lq . The theoretical levels are illustrated in each spin. The calculated E2 transition 
strengths are listed in Table IIIII along with experimental data [2^ . 

The calculated low-lying states have large overlaps with the basis wave functions in the small deformation region. 
For instance, the 3j" state has 87% overlap with (/3 cos 7,/? sin 7) = (0.20,0.09) (Fig. HJb)), which is the minimum 
point in the 3~ energy surface. As mentioned before, this state is interpreted as the intermediate between the 
cluster structure and the shell-model structure with the almost p3/2-subshell closed-proton configuration and the Ip- 
Ih excited neutron configuration. The calculated strength, B{E2, — )- Ij") = 1.3 e^fm^ is reasonable compared with 
the experimental value, B{E2,'i^ — > 1^) = 3.0 ± 1.2 e^fm**, though the level ordering is somehow in disagreement 
with the experimental one. As seen in the E2 transitions in Table IIIII E2 transitions among the low-lying states 
below the ^^He+a threshold are not strong and they show no remarkable collectivity nor specific band assignments 
for these states. 

In the highly excited negative-parity states above the ^°Be-|-a threshold energy, we obtain many developed cluster 
states having significant overlaps with the basis wave functions in the large /? regions. Some strong E2 transitions 
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TABLE I: Electromagnetic transition strengths B{E2) which are stronger than 1.0 e^fm"* for the positive-parity states in ^''C. 
The unit is e^fm*. 



Transition - 



Strength 



Theory Experiment 



2^ 
2^ 
2j 



2+ 

^5 



2+ 
2+ 



2| 
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2+ 
2+ 
2+ 



3^ 
3^ 
3^ 
3J 
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0^ 

ot 
ot 
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2+ 
2+ 

23+ 
2+ 



2J 

■2^ 
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67.2 
5.3 
1.6 
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2.8 
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2.6 
1.3 
6.3 
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2.5 
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1.7 
3.2 
2.0 



3.6 ±0.6 
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4^ 
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4^ 
4J 
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44+ 



2+ 
2+ 
2+ 

^3 



2+ 
2+ 
2+ 
2+ 

^8 

2+ 

^2 

2+ 

^4 

2+ 
2+ 

^8 

2+ 

^4 

2+ 
2+ 

^7 



■3+ 
■3+ 
■3^ 
3j 

■4| 
4 



■43+ 
■44+ 



4.7 
2.8 
2.1 
1.0 
5.1 
2.4 
12.9 
1.4 
2.5 
1.4 
4.4 
5.8 
1.0 
88.3 
8.7 

1.2 
4.6 
3.5 
2.8 

4.5 
1.7 
1.5 
3.0 
6.0 
2.1 



5+ ^ 3+ 2.4 



5f 
5^ 



■4+ 
■it 



4.9 
4.5 



are found in the calculated B{E2) values for high-lying states. However, the E2 strengths are fragmented into many 
states; therefore, it is difficult to classify these states in simple band assignments. In Fig. [3 we show the calculated 
level scheme of the high-lying negative-parity states above the ^"Be-f a threshold energy with E2 transition strengths. 
Among these highly excited states, some flows of strong E2 transition strengths are found in the present calculations. 
Transitions from each state are, however, dispersed. For example, the 5^ state has strong transitions to the 3^, 3^q, 
and 3j~2 states. These states are constituted by largely deformed states such as those shown in Figs. UJ^e), (d), and 
(e). What is different in the linear-chain band in the positive-parity states is that the GCM amplitudes of these 
negative-parity states do not concentrate on a single basis wave function with a specific geometric configuration but 
distribute into many basis wave functions. These features indicate strong mixing of the largely deformed prolate, 
oblate, and triaxial states. This mixing may be an origin of the complicated level scheme of the negative-parity states. 
As a result, the parity partner of the linear-chain band suggested in the positive-parity states may disappear in the 
negative-parity states. In spite of the largely scattered E2 strengths, there are sets of states classified by flows of 
significant strong E2 transitions. For example, a set consisting of 5^, 3j~Q, 3j~2, 3^3, 1^, 1^, and 1^ states can be 
interpreted as members of a quasiband with ^"^Be+a cluster structure in the negative-parity states of ^^C. In the next 
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TABLE II: Root-mean-square radii for mass distributions of the positive-parity states in ^ C. Tlie unit is fm. 
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FIG. 6: Energy levels of the negative-parity states in ^^C. Five columns on the left are the experimental data and five columns 
on the right are the calculated results. The dotted lines on the left and right show the experimental and theoretical ^"^Be-I-Q 
threshold energies, respectively. 



section, we explain a possible reason why this mixing occurs. 

We here comment the possible reason for missing of the low- lying 0~ state in the present calculation. It is expected 
to have Ip-lh configuration with one neutron in the lsi/2 orbital. In the present framework, the width parameter is 
taken to be a common value for all the single-particle Gaussian wave functions. The fixed Gaussian wave packets may 
not be suitable to describe the lsi/2 orbital because the spatial extent of the lsi/2 orbital should be large compared 
with other orbitals. 



IV. DISCUSSION 



In this section, we discuss the linear-chain band and the role of the excess neutrons. We also discuss the cluster 
structures of ^'^C in comparison with the findings of earlier works. 

First, wc discuss the reasons why the linear-chain band appears in the positive-parity states of ^"^C. For appearance 
of the linear-chain band, both the existence of the plateau around the linear-chain structure in the energy surface 
and the orthogonality to the lower states are essential. To explain more details, we show the squared overlaps 
|(0"'"|Po{j^°|$"'"(/3, 7))P of the GCM wave functions for the 0'*" states with the basis AMD wave function at each point 
(/3, 7) in Fig. [51 As mentioned before, the flat region exists near the largely deformed prolate region in the energy 
surfaces, and there exists a local minimum. Although this is a shallow local minimum with 1 MeV depth in the 0^ 
energy surface, the GCM amplitudes for the linear-chain 0^ state localize around this local minimum as shown in 
Fig. [H Another reason is the orthogonality to the triaxial bands consisting of the K"^ = J band and its K"^ = 2+ side 
band. In the case of 0+ states, the GCM amplitudes for the state occupy the triaxial region. We remind the reader 
that the bending 3a configurations are contained in the triaxial region. Therefore, even though the energy curve 
around the linear-chain structure is soft against the bending mode, the GCM amplitudes for the linear-chain 0'^ state 
are confined in the prolate region to satisfy the orthogonality to the triaxial state resulting in the concentration of 



12 



TABLE III: Electromagnetic transition strengths B{E2) which are stronger than 1.0 e^im^ for the negative-parity states in 
"C. The unit is e^fm"*. 
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the amplitude at the linear-chain structure. In other words, the triaxial state prevents the linear-chain 0^ state 
from bending. The situation is almost the same in the 2+ and 4+ states except that there are non-zero K components. 
In the 2+ case, distributions of the squared overlap on the /3-7 plane for the 2^ and 2^ states are similar to the 0^ 
and 0^ states, respectively. Moreover, the distribution of the squared overlap for the 2j^ state is similar to those for 
the 2^ state. This means the 2^ and 2^ states already exhaust the K'^ = 0+ and 2+ components of the bases which 
have bending configuration. Therefore, the 2+ states can be discussed in the same line as the 0''' states. That is, in 
the stabilization of the linear-chain state, the orthogonality of the 2^ state to the 2^ and 2^ states plays the same 
role as that of the 0^ state to the state. The 2j and 2^ states prevent the linear-chain 2^ state from bending, so 
that the GCM amplitudes for the 2^ state concentrate around the linear-chain structure. 

Next, we consider the role of the excess neutrons in the linear-chain structure of ^^C. The linear-chain Sa structure 
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FIG. 7: (Color online) Calculated energy levels of the negative-parity states above the ^"Be-I-Q: threshold energy with E2 
transition strengths in ^''C. An integer i at the right of a level denotes the order of states. 5.0 e^fm* < B{E2) < 10.0 
e^fm*, 10.0 e^fm'' < B{E2) < 20.0 e^fm", 20.0 e^fm" < B{E2) < 50.0 e^fm*, and 50.0 e^fm* < B{E2) transitions are described 
by broken, black solid, green bold solid, and red bold solid arrows, respectively. The dotted line shows theoretical ^^Be+a 
threshold energy. 



in C isotopes has been a long standing problem, for example, as discussed in Refs. [3, [l^, [U [s^l ■ In the recent works 
on ^^C with AMD and FMD calculations [1, [T^jH^Ij open triangle 3a structure was suggested only in the Oj^ state, 
but no straight-line 3Q;-chain structure nor linear-chain band was found in the ^^C system. It is contrastive to ^^C 
where the straight linear-chain structure is suggested in the 0^, 2^, 4^ states in the present work. It means that 
there is no linear-chain state in ^^C except for the 0;^ state, while the linear-chain structure construct the rotational 
band in ^^C. By comparing the present results of ^^C with the ^^C results calculated with the same framework of 
the /3-7 constraint AMD -I- GCM in Ref. [l^l, it is found that the difference between the and the ^^C cases 
originates in the differences of the behavior of the energy surfaces and also the structure of the states lower than 
the linear-chain structure. First, we describe the difference in the energy surfaces. In Fig. 1^1 we show the energy 
surfaces for ^^C to compare them with the those for ^''^C. The detailed discussions for ^^C are written in Ref. [20l |. 
In the energy surfaces of ^^C (Fig. 9), there is a flat region around the largely deformed prolate region but there is 
no local minimum even in the 0"*" energy surface. However, in the energy surfaces of ^^C (Fig. 1). there is a shallow 
local minimum with 1 MeV depth around the largely deformed prolate region. This difference indicates that excess 
neutrons give the effect to favor the linear-chain structure. Second, we discuss the structure of the lower states than 
the linear-chain structure. As explained in the previous paragraph, existence of the triaxial bands is essential for 
the linear-chain band in ^^C because they prevent the linear-chain band from bending because of the orthogonality. 
However, in ^^C, the states which prevent the linear-chain states from bending do not exist. An exception is the case 
of 0+ states. Let us first explain the 0+ cases, which are similar in ^^C and ^*C. In ^^C, GCM amplitudes of the 
state occupy the considerably bending 3a structure and distribute widely in various basis wave functions. Therefore, 
the 0^ state prevents the 0^ state from bending extremely because of the orthogonality. Because the lower state, Oj, 
already exhausts the bases which have bending configuration, the higher state, 0^ , cannot have bending configuration; 
therefore, the linear-chainlike structure dominates the 0'^ state even though the linear-chain structure is not a local 
minimum in the energy surface in ^^C. This is similar to the relation between the 0^ and 0^ states in ^"^C. However, 
for the 2+ and 4+ states in ^^C, there is no such a state that exhausts bases with bending configurations. Therefore, 
mixing of the largely deformed prolate, oblate and triaxial states occurs. To summarize, the excess neutrons play 
two roles in the appearance of the rotational band of the linear-chain structure. First, because of the existence of the 
excess neutrons, the energy surfaces for the positive-parity states change to favor the linear-chain structure. Second, 



14 




1.0 
0.9 
H0.8 
0.7 
H0.6 
0.5 
0.4 
0.3 
0.2 
0.1 

o:o 



0.6 0,3 



FIG. 8: Overlaps of the GCM wave functions for the 0^ states of ^''C with basis AMD wave functions, Pqq 7)), obtained 

with the /3-7 constraint. The overlap at a certain {{3, 7) point is shown in the contour map. 



the excess neutrons are essential for constructing the triaxial bands, vifhich occupy the bending 3a configurations 
considerably and prevent the linear-chain states from bending to satisfy the orthogonality to each other. 

As we explain in the previous section, the linear-chain band does not appear in the negative-parity states. The 
reason is clear from the behavior of the energy surfaces for the negative-parity states. There is no flat region in 
the largely deformed prolate region. This means the ^"Be cluster rotates easily in the negative-parity states. This 
nature was already explained also in the discussion of the intrinsic density distributions shown in Fig. SJ^e). Because 
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FIG. 9: Energy surfaces of C on the /?-7 plane. The upper shows the energy for the positive-parity states before the 
total-angular-momentum projection and the lower shows that for the 0+ states after the total-angular-momentum projection. 



of this unstable feature against rotational mode of the Be cluster, the prolate deformed states strongly mix with 
the triaxially deformed states in the negative-parity states. 

Linear-chain structures of C isotopes have been discussed in earlier works with 3a+Xn cluster models [HI, [l^, [l^ . 
In the study by Itagaki et al. [l8| . it was reported that the linear-chain state may not be stable in ^*C. Their argument 
seems to contradict the present results of the linear-chain band suggested in the positive-parity states of ^"'C. In their 
work, it was shown that the linear-chain structure of ^"'C is unstable against the bending mode because there is no 
local minimum in the energy curve. However, to confirm the stability of an excited state, it is essential to take into 
account the orthogonality to lower states as done in the GCM calculations. In the results of Ref. [l^, the energy 
curve along the bending mode is very flat and is not so much different from the present results though it has no 
local minimum in their calculations while the shallow local minimum exists in the flat region of the present results. 
As mentioned before, the linear-chain band appears from this flat region mainly because of the orthogonality to the 
triaxial bands which block the bending of the 3a configuration. Therefore, we should stress again the importance of 
the GCM calculations for excited states. Another difference of the model used in Ref. [l^ from the preset framework 
is that three a clusters are set at the same intervals and the excess neutrons are assumed to occupy the orbital 
moving around the whole 3a core. In such a model, the system has to have a parity symmetric structure and the 
^"Be correlation cannot be incorporated. In the present calculations, though the 3a core is not a priori assumed in 
the framework, various 3a-|-2n structures are obtained in the energy variation. The geometric configuration of three 
a clusters and Gaussian centers for excess neutrons are optimized. As a result, the linear-chain states show the ^"Be 
correlation as discussed before. In further detailed analysis of the linear-chain structure, we confirm that the ^"Be 
correlation gains energy compared with the symmetric configuration of the linear-chain 3a-t-2n structure. These will 
be reported in a future article. 

The idea of ^°Be-|-a structure in ^"^C has been proposed by von Oertzen et al. in Ref. [l2|. They conjectured that the 
parity doublet band could be built from the parity-asymmetric linear-chain structure because of the ^"Be correlation. 
However, the present results suggest that the linear-chain structure gives only the positive-parity rotational band 
but no negative-parity band. As described before, because the ^"Be core can easily rotate in the negative-parity 
^°Be-|-a states, the strong mixing of the linear-chain structure with various bending -'^''Be-l-a configurations occurs; 
therefore, the parity partner of the positive-parity linear-chain band disappears. Another difference between the 
present results and their suggestion is the energy position of the positive-parity linear-chain band. The linear-chain 
band is obtained above the ^°Be-|-a threshold in the present calculations, whereas they proposed the linear-chain 
band below the threshold. The reason for the higher energy than the ^"^Be-l-a threshold can be naturally understood 
by kinetic energy loss to align ^"^Be and a clusters in the linear-chain configuration. 

Recently, Itagaki et al. have studied excited states of ^"'C with a 3a-t-2n cluster model and have predicted the 
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equilateral-triangular structure of the developed three a clusters with excess neutrons [17[ . They have proposed that 
the excess neutrons moving in molecular orbitals around the three a clusters stabilize the geometric 3a configuration 
of equilateral triangle and suggested two rotational bands, K'" ~ O2 and K'" — 3~ , resulting from the equilateral 
triangle intrinsic structure. For the positive-parity states, we also obtained the = band with well-developed 
three a-cluster cores in an almost equilateral-triangular configuration, and therefore, in a sense, the present results 
seem to be partially consistent with the work of Itagaki et ai. However, the present results suggest the triaxiality 
owing to the excess neutron orbitals and the resultant K"^ = 2+ side band on the K"^ = band. These are 
considerably different from the earlier work. In addition, the development of three a clusters is weaker than that 
suggested in Ref. [l3| as the root-mean-square radius for the state is 2.45 fm in the present calculation while it 
is about 2.6 fm in the earlier work. For the negative-parity states, a K"^ = 3~ rotational band starting from the 
3^ state was suggested in Ref. However, any clear = 3~ rotational band cannot be found in the present 

calculations. As mentioned, our results suggest that the a breaking is rather important in the low-lying states of ^''^C. 
The description of the a breaking is one of the advantages of the present model, while molecular orbitals for excess 
neutrons incorporated by Itagaki et al. are not completely included in our model. It is requested to investigate 
in more details with a model including both of the a breaking and molecular orbitals around a clustering core. 



V. SUMMARY AND OUTLOOK 



We investigated structures of the ground and excited states of ^"^C with the method of /3-j constraint AMD + 
GCM. The results reproduce well observed data, such as the energy levels except for 1+ and 0^ states, the root-mean- 
square radii, and the E2 transition strengths. The present results suggest that the ground state has the intermediate 
structure between cluster and shell-model structures. In the excited states, well-developed cluster structures are 
found in both the positive- and negative-parity states. Owing to the 3a configurations and the excess neutron motion, 
various structures appear. 

By analyzing the E2 transition strengths as well as the GCM amplitudes, we assigned members of four positive- 
parity bands: the ground, the triaxial = Oj, the triaxial = 2+, and the linear-chain bands. The K"^ = band 
and the K"^ = 2+ side band are constructed from the triaxially deformed intrinsic state, in which three a clusters have 
an almost equilateral-triangular configuration and the excess neutrons occupy the sd-like orbital. In the linear-chain 
band, the excess neutrons distribute around two of the three a clusters, which indicates the ^"^Be-l-a correlation. This 
is the first work that suggests the triaxial bands in ^'^C. For the 3a linear-chain band, which has been an attractive 
topic in C isotopes, there have been few theoretical works suggesting the appearance of the linear-chain band in ^*C. 
In the present work, it was found that excess neutrons play two important roles in the appearance of the rotational 
band of the linear-chain structure as follows. First, because of the existence of the excess neutrons, the energy surfaces 
for the positive-parity states change to favor the linear-chain structure. Second, the excess neutrons are essential for 
constructing the triaxial bands, which occupy the bending 3a configurations considerably and prevent the linear-chain 
states from bending to satisfy the orthogonality to each other. This result suggests an interesting mechanism for the 
stabilization of the exotic cluster structures owing to the excess neutron effects. 

For the negative-parity states, because of the largely scattered E2 strengths, we could not assign simple band 
structures. The origin of these dispersion of the E2 strengths is the strong mixing of the largely deformed prolate, 
oblate, and triaxial states in the negative-parity states. Although the positive-parity linear-chain band with the 
"'^'^Be-l-a correlation was suggested, its parity partner with negative parity disappears because of this mixing. In spite 
of the largely scattered E2 strengths in the negative-parity states, there are sets of states with significant strong 
E2 transitions. These states can be interpreted as members of a quasiband with ^"^Be-l-a cluster structure in the 
negative-parity states of ^''C. 

As mentioned earlier, the linear-chain structure shows the ^''Be correlation with the asymmetric structure instead 
of the symmetric configuration of the linear-chain 3a-|-2n structure. Further analysis of the ^"Be correlation in the 
linear-chain structure will be reported in a future article. 

The theoretically suggested excited states with ^°Be-|-a cluster structures are the candidates for the excited states 
which were observed recently in ^^Bc+a decays. The decay widths and the branching ratios to ^°Be(0^)4-a and 
^'^'Be{2f)+a could be helpful for assigning the calculated states to the observed states. Theoretical estimation of the 
partial decay widths of the excited states above the ^"Be-l-a threshold is also a remaining future problem. 

In the present calculations of ^^C, although a clusters are not assumed in the framework, three a-cluster cores 
are found to develop in excited states. Owing to a variety of 3a configurations and excess neutron motion, various 
interesting structures appear in the "'^^C system, where the excess neutrons play important roles in the appearance of 
cluster structures. It is a challenging problem to investigate the excited states of further neutron-rich C isotopes such 
as ^^C while focusing on the possibility of cluster states with three a-cluster cores with more excess neutrons. 
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